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ABSTRACT 

We use a new method to trace backwards the star formation history of the 
Milky Way disk, using a sample of M dwarfs in the solar neighbourhood which 
is representative for the entire solar circle. M stars are used because they show 
Hq, emission until a particular age which is a well calibrated function of their 
absolute magnitudes. This allows us to reconstruct the rate at which disk stars 
have been born over about half the disk's lifetime. Our star formation rate agrees 
well with those obtained by using other, independent, methods and seems to rule 
out a constant star formation rate. 

The principal result of this study is to show that a relation of the Schmidt- 
Kennicut type (which relates the star formation rate to the interstellar gas content 
of galaxy disks) has pertained in the Milky Way disk during the last 5 Gyr. The 
star formation rate we derive from the M dwarfs and the interstellar gas content 
of the disk can be inferred as a function of time from a model of the chemical 
enrichment of the disk, which is well constrained by the observations indicating 
that the metallicity of the Galactic disk has remained nearly constant over the 
timescales involved. We demonstrate that the star formation rate and gas surface 
densities over the last 5 Gyrs can be accurately described by a Schmidt-Kennicutt 
law with an index of F = 1.45^9 09. This is, within statistical uncertainties, the 
same value found for other galaxies. 

Subject headings: Galaxy: formation, Galaxy: evolution 
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1. Introduction 

Tracing the star formation history (SFH) in galactic disks is a crucial aspect of 
understanding their evolution. Determination of a detailed SFH is practically only possible 
in the disk of the Milky Way, and va rious me t hods have been developed to do so. The 
most common method, introduced by iTwarog] (119801 ). is to model the colour - magnitude 



diagram (CMD) of stars in the solar neighbourhood. Model CMDs are computed by 
adopting an initial mass function, a metallicity - age relation, sets of isochrones and 



iterating the S 



BerteUi fc Nasi 



H until computed CMD s match obseryation s ( Hernandez et al. 



2001 



Vergely et al 



2002 



Cignoni et al. 



20061). 



200G 



Rocha-Pinto et al 



(l2000h 



used the chromospheric emission intensity in the Ca H and K absorption lines in the spectra 
of F and G stars. However, due to stellar activity cycles, metallicity effects and other 



complications, chromospheric emission dating can become 



2001 



for a full discussion). A third method, described by 



very uncertain (see 



Noh fc Scalol Jl990h. is based 



Feltzing et al. 



on 



the lu minosity function of white dwarfs, which preserves details of the SFH. iLada fc Lada 



( 120031 ) have demonstrated that a very considerable fraction of the field star population in 



the Milky Way disk might have been formed in open clusters. These dissolve on fairly s. 



lort 



timescales due to dynamical effects ( jde la Fuente Marcos fc de la Fuente Marcos 



2004 ) 



so that, even though only a small percentage of the stellar population is actually found 
in open clusters, the SFH of open clusters might be indicative for the overall SFH of the 
Galaxy. However, since open clusters are typically young, this method can be used to study 
the SFH only over the last Gyr of the evolution of the Galactic disk. Despite differences 
in the details, the broad result of these studies is that the star formation rate (SFR) has 
been declining slowly for the last 5 Gyr, prior to which it was rising slowly, and has been 
punctuated by bursts of star formation from time to time. 



In this study we analyze a sample of M dwarfs, whose spectra show Hq, emission 
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(dMe), which we have drawn from the Fourth Edition of the Catalogue of Nearby Stars 



(CNS4, 



JahreiB 



20081 ). Like chromospheric Ca H and K emission, Hq, emission is related to 



magnetic activity in the stars which is in turn regul ated by their inte r nal constitution, and 



2004). 



Hawley et al. 



(1200 



in particular by the depth of their convection zones (IWest et al. 
have demonstrated, using M dwarfs in nearby open clusters, that the period in which the 
stars show emission is (among other stellar properties) closely correlated with their 
absolute magnitudes, and ranges from about lO"^'^ yrs at absolute \^-band magnitude My = 
6 to some 10^'^ yrs at an absolute magnitude of My = 13. After these timescales, the Hq, 
emission essentially switches off. We have utilised this effect to reconstruct the SFH of the 
Galactic disk in a new, independent manner. 

In galaxy evolution studies, the SFR is commonly expressed as a "Schmidt-Kennicutt" 
power law function of the surface density of the interstellar gas, e ven though the physical 



proce sses leading to this relation are not yet fully understood (cf. 



2008 



Schaye. fc Dalla Vecchia 



for a discussion of the various concepts). A relation betwe en the S 



Schmidt 



-j and the density 



(119591 ) in order 



of interstellar gas raised to some power was first suggested by 
to explain that the vertical scale height of young stellar populations in the Galactic disk 
is smaller c ompared to that of the interstellar gas. This relation was later modified by 



Kennicutt 



(11989 



19981 ) who correlated the SFR observed in other spiral galaxies with the 
surface densities of their interstellar gas either in concentric rings around the centers of the 
galaxies or averaged ov er their entire optical d isks. The SFRs and gas surface densities of 



the galax ies studied by 



Similarly iBoissier et al. 



Kennicutt 



(11989 



19981 ) are based on present day states of the disks. 



(I2OOII ) have derived a Schmidt-Kennicutt relation for the disk of 
the Milky Way based on its present day gas content. They correlated the radial density 
distribution of the interstellar gas in the Galactic disk with the radial variation of the disk's 
present day SFR, finding a relation which is steeper than the original Schmidt-Kennicutt 
relation for other galaxies. What we address here is whether the Schmidt-Kennicutt law 
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also holds true in the disk's evolution in the past, particularly over the last 5 Gyrs or so. 
The M dwarfs permit us to track the SFH, and we use the chemical enrichment history of 
the Galactic disk to reconstruct its gas surface density over time. 

Our paper is organized as follows: In the next section we describe our data. In section 
3 we derive the SFH and in section 4 we try to establish a Schmidt-Kennicutt law for earlier 
epochs of Galactic evolution. The conclusions are summarized in the last section. 



2. Data 



Our basic sample consists of 1453 M dwarfs extracted from the CNS4 (4th Catalog 
of Nearby Stars) which are within 25 pc from the Sun. 210 of them are classified as 
dMe. T he emission strength data come frona the Palomar/MSU Nearby-Star Spectroscopic 



Survey (IReid et al. 



1995 



Hawlev et al 



of the Catalogue of Nearby Stars (iGliese fc Jahreifi 



1996r). This survey was based on an earlier version 



19911 ). The dMe are binned into 
half-magnitude intervals of absolute magnitudes from My = 8.75 to My = 12.75 mag and 
are shown in Fig. [1] together with star counts in the CNS4 of all M dwarfs of the same 
absolute magnitudes. The latter were restricted to that sample on which the Palomar/MSU 
Nearby-Star Spectroscopic Survey is based. Since we want to concentrate here on the thin 
disk of the Milky Way, we only use stars with vertical W velocities with \W\ < 50 km/s 



( iNordstrom et al 



2004h 



To assess the completeness of our sample, we have determined the cumulative number 
of stars as a function of distance from the Sun. These are shown in Fig. [2] on double 
logarithmic scales. In a spatially homogeneous sample the cumulative number of stars grows 
in proportion to r^. This is indicated as solid lines with a slope of 3 in the logarithmic 
diagram. As can be seen from Fig. [2] our samples become seriously incomplete beyond 15 
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pc of the Sun (this was demonstrated aheady by lJahreiB fc WielenI 119971 with a prehminary 
version of the CNS4), but this incompleteness has only a small effect on the method we 
use, which depends on the relative numbers of stars with and without H„ emission, and 
furthermore such stars have almost the same likelihood of being in the source catalogue. 
This can be quantified as follows. Fig. [2] indicates that the dMe are undersampled compared 
to the entire sample of M dwarfs. For example, within a radius of 25 pc from the Sun 53 
percent of nearby M dwarfs in the magnitude range, which we consider here, are catalogued 
in the CNS4, and 44 percent of the dMe in this range are so catalogued. To reach a 
completeness level the same as in the complete sample at this distance (25 pc) implies that 
some 43 dMe are missing from the CNS4. 

The different sampling rates are almost certainly due to the fact that the M dwarf 
sample in the CNS4 is biased towards high proper motion stars. dMe at the bright end 
of the luminosity function are very young and have low space velocities. Thus it is very 



plausible that 



Hawley et al. 



they are undersampled compared to the entire M dwarf sample. Indeed 



(119961 ) have shown that the velocity dispersions of the dMe are significantly 
smaller than the mean velocity dispersions of M dwarfs in the CNS4. We discuss possible 
consequences of the effect for our results below. 

We note that dMe seem frequently to be members of binary or multiple systems. 63% 
of the dMe of our sample are flagged as binaries in the CNS4, whereas only 28% of the M 
dwarfs without emission are so flagged. However we do not attribute any significance 
to these numbers, because the degree of co mpleteness of the identifi cation of binaries in the 



CNS4 is not known. In a systematic study. 



Fischer fc Marcy 



(11992h found that 42%±9% of 



M dwarfs in the solar neighbourhood are binaries. Unresolved binarity might occasionally 
shift a star in our sample from one half-magnitude bin into a neighbouring bin, but we 
expect no systematic effects as a consequence. 
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Next we calculate the li mitin g age Tumit of the stars in each magnitude bin, using the 
calibration of 



Hawley et al. 



(120001) 



-16.75 + 3.06 logio(ri. 



imit ) 



(1) 



Star formation rates are usually reckoned in terms of surface density per time interval, 
whereas our sample is confined to within 25 pc of the Sun, and samples volume density. 
To correct our sample from volume to surface density, we note that in the disk's vertical 
gravitational potential the vertical scale height of stars is known to depend linearly on their 
vertical velocity dispersion. We have multiplied the numbers of M dwarfs in each magnitude 
bin by the W velocity dispersion of stars with the same typical ages to obtain the surface 
density of stars of that age, N\Mv) = N{Mv) ■ awi^). For the latter we have adopted 



the va lues given for groups 6d, 6c, 6b, 6a, and groups 5 - 2 in Table 4 of 



JahreiB fc Wielen 



fll997h . 



Our derivation of the SFH relies thus on the frequency of M dwarfs observed in the 
solar neighbourhood and corrected to be representative for a cylinder at the position of the 
Sun and perpendicular to the Galactic plane. However, the local volume actually samples 
stars from much more distant parts of the Milky Way. At the solar annulus, the orbital 
period around the Galactic center, p = 27t/Qq, where Qq denotes the angular velocity of the 
local standard of rest, is p = 2tt/ 220 km/s / 8 kpc = 2.3-10^ yrs. This means that stars 
with ages of order 10* yrs can reach the solar neighbourhood from the other side of the 
Galaxy. This effect is enhanced by the epicyclic motions of the stars. The typical sizes of 
the epicycles can be estimated from the velocity dispersions of the stars txr/ in th e radial 



and (Ty in the azimuthal directions, respectively, with the relations (jWielen 



19821) 
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c^K=t/r|tB+^ and cr^ = ^--R. (2) 
B denotes the second Oort constant. If we assume a flat Galacti c rotation curve, 



B = — K^n, and use again the velocity dispersions given in Table 4 of 



Jahreifj Sz Wielen 



(119971 ). we find for stars with ages r = 0.2 Gyrs = 240 pc and cr$ = 340 pc, rising to aji 
= 1.1 kpc and cr$ = 1.6 kpc for stars with a mean age of r = 2.3 Gyrs, which correspond to 
our oldest age bin. Thus we sample dMe in a ring around the Galactic center with a width 
of about one to two kpc, which we term "the solar circle" . 



3. Star formation history 

The surface densities determined in the previous section depend on how far back 
the past SFH has been sampled, but are also weighted by the Main Sequence luminosity 
function ^rnsiMy), i.e. the relative frequency distribution of stars as a function of absolute 
magnitude. Since M dwarfs are so long-lived, their luminosity function can be safely 
assumed to be unchanging with time. As is customary practice in studies of the SFH we 
assume that the SFH has been the same for all M dwarfs in our s ample irrespective of mass 



(they have masses in the range 0.7 to 0.2 Mq, iHenry et al.lll993l ). Below we compare our 
results with SFHs determined with other stellar types. The numbers N'{Mv) determined 
above can be then interpreted as 



OC 



S.(T) ' 

where T is the present age of the Galactic disk, for which we adopt a value of 10^" 



(3) 
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yrs, and where and denote the surface density of the stellar disk and the SFR 

at the epoch t = T — TumuiMy), respectively. In eq. Owe have assumed that the surface 
density of stellar disk was initially S^(0) = 0. Thus the ratio N'^^{My)/N^^^{Mv) allows 
us to determine S*(T) — which represents the SFH of the disk as a function of 

time t relative to the stellar surface density now {t = T). This is shown on logarithmic 
scale in Fig. [31 The vertical error bars reflect only Poisson errors. The horizontal error 
bars have been estimated by th e scatter of th e calib ration data points around the linear 



My - log (age) regression line of 



Hawley et al 



(I2OOOI ). The calibration is based on 5 open 
clusters in which the authors have determined the limiting My of the dMes. The ag es of the 



clusters are somewhat uncertain. To their oldest cluster, M67, 



Hawley et al 



an age of 5.4 Gyrs, whereas ot her determinations find 4 ± 0.5 Gyrs ( 



(120001) 



Dinescu et al 



Mamaiek fc Hillenbrand 



20081 ). The ages of the younger clusters used by 



assig n 



1995 



Hawley et al 



( I2OOOI ) like the Hyades or the Pleia desof with ages of 625 and 130 Myrs, respectively, are 



uncertain on the order of 10 Myrs (IMamajek fc Hillenbrand 



20081). 



We have fitted to the data in Fig. [3] in the range of disk age t = 5 to 10 Gyrs a relation 
of the form 



(4) 



and find an index of = 2.2_Qg. The star formation rate is then given by 



-(1+/3) 



(5) 



a relation which declines with time. As can be seen from Fig. 3 a constant star 
formation rate does not fit the data. 



As described above, the ratios N'^ /A^tot have to be corrected upwards due to the 
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relative undersampling of the dMe. This affects mainly the younger, bright dMe as they 
have low space velocities. In our view this explains the irregular pattern of N'^^/N[^^ ratios 
at t = 9 and 9.3 Gyrs. Shifting the latter up to the fitting curve in Fig. [3] would require 
some 39 additional stars. This is almost precisely the number of stars we have estimated 
above needed to be added to the dMe sample so that it reaches the same level of relative 
completeness as the entire M dwarf sample. As can be seen from Fig. [3] such a correction 
leads to an even better fit of eq. IHto the N'^^/N[^^ data. 

The Sun and the stars in its vicinity are at present embedded in the "Local Bubble" 
of the interstellar medium, a cavity with a diameter of about 200 pc filled with hot diluted 
gas. The physical conditions of this gas are such that one does not exp ect any stars to be 



Fuchsetal 



(120061) 



presently forming within it. However, the Local Bubble is quite young, 
have recently redetermined its age and found that it was formed about 15 Myrs or less ago. 
The ages of the dMe in our youngest age bin are of the order of 100 Myrs so that most of 
them were born before the Local Bubble was formed. Our results on the SFH are therefore 
not affected by the presence of the Local Bubble today. 

Unfortunately we cannot derive absolute rates for the SFH using our method. However 
we can estimate its magnitude roughly by observing that, according to our determination, 
the surface density of the stellar component of t he Galactic d i sk's m ass grew by 20 percent 



of today's surface density over the last 3 Gyrs. 



Flynn et al. 



(120061 ) find a local surface 



density of the stellar Galactic disk of 35.5 M.q pc ^ which corresponds to a SFR of 2.4 



Mq pc-2 Gyr-1. The 
is 6.5 M.Q pc~^ ([Dame 



ocal s urface density of interstellar atomic and molecular hydrogen 



1993r). Very in t eresti ngly, this combination of parameters places 



the Milky Way right on the iKennicutt 



(119981 ) relation, even though we use local disk star 



data, whereas Kennicu tt's data are radi ally averaged over the whole disk. In Fig. S] we have 



reproduced the data of 



Kennicutt 



(119981 ) and indicated the location of the Milky Way in the 
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diagram. As can be seen from Fig. |H the Milky Way falls precisely in the region populated 
by similar galaxy types to it, i.e. Sb-Sbc. 

Our determination of the shape of the SFH is consistent with other determinations of 



Cienoni et al. 


( 


2006[). 


de la Fuente Marcos & de la F 


uente Marcos (2004). 


Hernandez et al. 


(2000 




Rocha-Pinto et al. 


(2000 


), and 
of star 


Vereelv et al. 


(2002 


). The results of 


Hernandez et al. 


(2000 


who modelled the CMD 


3 vouneer than 3 Gvrs. and Vereelv et al. 


(2002 


), who 



modelled the CMD of nearby stars observed with Hipparcos, can be directly compared with 
our results, because they refer to column densities. In Fig. O we show on a logarithmic scale 



S^,(T) — as obtained in this work versus (T) — E 



by integrating the SFRs determined by 



Hernandez et al. 



(t), w hich we have cons t ructed 



(120001) and 



Vergely et al 



(120021) 



Short duration star bursts in the SFR are effectively smoothed away by the integration 
procedure. We show in Fig. 5 also a SFH de t ermin ed by integrating the SFR derived by 



de la Fuente Marcos &: de la Fuente Marcos 



(I2OO4J ) from the age distribution of nearby 



open clusters. Their local number densities have been corrected to be representative for 
column densities by multiplying them by the vertical velocity dispersions cr(r) typical for 



stars of the same age (IJahreifi fc Wielen 



the SFR found by 



Cignoni et al 



19971 ). Similarly we have integrated and converted 



( I2OO6I ) modelling the CMD of Hipparcos stars (using the 
data which have been clipped in velocity space). Over the last 3 Gyrs all SFHs fit excellently 
together. Given the totally different approaches in deriving the SFH the good agreement 
is, in our view, quite remarkable. As explained above, the deviation of our data 0.7 and 1 
Gyr ago from the other SFHs and the fitting formula H] is in our interpretation due to the 
proper motion bias in the CNS4. This affects mainly the bright dMes with short limiting 
in emission ages. The v elocity dispersion of young stars increases with age by factor of 3 



JahreiB &: Wielen 



19971 ). Thus our older age bins are expected to b e more comple t e as t hey 

( 20061 ) is 



indeed seem to be. Beyond ages of 3 Gyrs the SFH determined by 



Cignoni et al 
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consistent with the SFH derived in this work, whereas the SFH by lVergely et aLl (120021 ) falls 
in that age range for unknown reasons systematically below the two oth er determinations. 
Final 



Rocha-Pinto et al. 



y we show in Fig. 5 also a SFH obtained by integrating the SFR of 
( I2OOOI ) which is based on nearby F and G stars dated with chromospheric emission ages. 
This SFH differs significantly from the other SFHs shown ther e and is consistent with a 
constant SFR (cf. Fig. [3]). With the exception of the results by iRocha-Pinto et al.l (I2OOOI ). 
the consistent picture of the disk's SFR is that S*(T) — S*(t) rises steeply as a function of 
look-back time T — t until about 5 Gyr. We cannot resolve with our data the SFH at earlier 
epochs of Gal a ctic e volution. However, the other determinations of the SFH summarized in 



Cignoni et al. 



(120061 ) indicate a turn over of the SFR to an apparently much reduced rate 



at early epochs. 



4. Schmidt-Kennicutt law 



As 



Kennicutt 



( 1l998l ) has pointed out, the data (which we have reproduced in Fig. H]) 
suggest a relation between the SFR and the surface density of the interstellar gas of the 
form 



(6) 



with an index F = 1.4 ± 0.15 and a coefficient k = (2.5±0.7)-10"'' Mg/yr/kpc^. 

Let us consider the gas content of the Galactic disk as being regulated by two processes 
— on one hand interstellar gas is converted into stars, but on the other hand the disk 
accrets primordial gas. 



(7) 
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The present accretion r ate of gas, as inferred from in 



is about 1 A4q pc ^ Gyr ^ (IWakker fc van Woerden 



ailing high velocity HI clouds 



199ll ). but the rate is unknown for 



earlier epochs of Galactic evolution. In order to constrain the amount of accreted gas we 
invoke constraints provided by the chemical enrichment history of the Galactic disk. In a 
one-zone model the enrichment or depletion of heavy elements is described by 



dt 



[ZT.,) = -Zt, + yt, 



(8) 



where Z denotes the mass fraction of heavy elements. In eq. [HI we employed the 
instantaneous recycling approximation charac terized by the yield y . The Geneva - 



Copenhagen Survey of nearby F and G stars (INordstrom et al 



20041 ) has shown that the 



metallicity of the Galactic disk has been, within statistical uncertainties, con stan t over most 



of its e volution ("flat"). In vie w of contradictory claims by iHaywoodl (120061 ) and 



Reid et al. 



mm 



Holmberg et al. 



(120071 ) have thoroughly reanalyzed the Geneva - Copenhagen survey 
and discussed meticulously all selection effects. They confirm their earlier conclusion that 
the age-metallicity relation is indeed flat. Assuming this to be the case, equations [3] and [H] 



imply then S 



accr 
9 



{y/Z) or 



(9) 



The ratio y/Z has a value of about yjZ = 0.7 (jPagell 119971 ). which is another way of 
stating that accretion of gas is important for the evolution of the Galactic disk. 

Taking the derivative of eq. [6] we find 



(10) 
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If we insert the parameterization of introduced in eq. HI eqs. O [H and [TO] lead to 



with 



log 



consti + - log (Sg - Sgo) 

i + /3r-i 

consta H ^ ^ log (Sg - Sgo) 



E,o = S,(T) - (l - I) ■ - S.(T)) . 



A comparison of the coefficients shows then that the indices are related by 



l + P 
/3 



(12) 



In the previous section we determined a value of /3 = 2.2lQg, which implies an index of 
the Schmidt-Kennicutt law of F = lAb'^ n m- T his fits, within statistical uncertainties, very 



well to the index observed by 



Kennicutt 



19981 ) for the present day SFR in other galaxies. 



Apparently the SFH of the Milky Way disk was regulated by the same Schmidt-Kennicutt 
law over the last 5 Gyrs. At earlier epochs the slope (3 was probably shallower, which would 
imply a larger index F for the Schmidt-Kennicutt law. 

In deriving eq. [TT]we have assumed an accretion rate of primordial gas that compensated 
exactly the metal enrichment due to star formation so that the metallicity stayed constant 



during the evolution of the Galactic disk, in line with the observations by 



Holmberg et al.l 



(120071 ). But many models of the chemical enrichment history of t he Galaxy predict a mild 



Colavitti et al. 



2008 



for a 



chemical enrichment of the Galactic disk over the last 5 Gyrs (see 
recent discussion). This is due to less accretion than we assumed. However, our principal 
result on the slope of the Schmidt-Kennicutt law does not depend on the absolute amount 
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of accreted gas. We note that a 'closed-box' model with E^'^'^^ = 0, which is known to lead 
to a relative distribution of metallicities in the Milky Way disk quite inconsistent with 



observations (jPagel 



19971 ). implies the same relation as we found in eq. [T21 Cosmological 



models of gas accretion onto galaxies and of the ensuing chemi cal enrichment history predict 



accretion rates which decline with time (IColavitti et al. 



20081), similar to the SFR which 



we found for the last 5 Gyrs. Thus a relation of the kind E^^"' ~ S^, will approximately 
hold and we argue that our result on the slope of the Schmidt-Kennicutt relation should be 
rather robust. 



5. Conclusions 

We have used a new method on a sample of M dwarfs in the solar neighbourhood 
showing Hq, emission in their spectra, to trace backwards the star formation history (SFH) 
of the Milky Way disk. The M stars show Hq, emission up until a certain age (which is a 
well calibrated function of their absolute magnitudes) allowing us to reconstruct the rate at 
which stars have been born in the disk at the solar circle over about half its lifetime. 

The star formation history of the Galaxy seems to be well resolved over the last 5 Gyrs 
using our M dwarf sample. We find that the star formation rate (SFR) rises continuously 
when looking backwards to 5 Gyr. This is in very good agreement with other determinations 
of the star formation history of the Galactic disk by different methods. 

In order to derive a Schmidt-Kennicutt law for the local disk, we have determined 
the variation of the surface density of the interstellar gas in the Galactic disk during the 
formation and evolution of the disk. The gas content of the disk is controlled by the 
competition between the star formation process, on one hand, and accretion of gas from 
outside the Galaxy, on the other hand. Since the latter cannot be observed directly, we had 
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to infer the amount of accreted gas indirectly. For this purpose we have invoked the chemical 
enrichment history of the Galactic disk. Based on the observation that the metallicity of 
the Galactic disk has remained rather constant during its evolution (in particular over the 
last 5 Gyr) we show that the gas infall rate is proportional to the SFR. 

We have then correlated the SFR with the surface density of the interstellar gas and 
found that the correlation followed over the last 5 Gyrs a Schmidt-Kennicutt law with 
an index of F = l^Sl^ T his ag rees, within statistical uncertainties, with the value of 



F = 1.4 found by 



Kennicutt 



(119981 ) for the present day SFR in other galaxies. It seems 
that at epochs of Galactic evolution earlier than that stars formed at a much reduced star 
formation rate, which does not follow the Schmidt-Kennicutt law. The investigation of the 
reason for this change of the mode of star formation is beyond the scope of this paper, but 
we hope to address this question in the future. 
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Fig. 1.— Star counts of M dwarfs in the Catalogue of Nearby Stars ( iGliese fc Jahreifi 



13 



1991 



JahreiB 



20081 ) as function of absolute magnitude. The light blue histogram shows the number 
of 210 M dwarfs with Hq, emission and the dark blue histogram shows the distribution of all 
1453 M dwarfs in the magnitude range My = 8.75 to 12.75 mag, respectively. 
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Fig. 2. — Radial cumulative number distribution of M dwarfs with absolute magnitudes My 
= 8 to 13 mag within the solar neighbourhood drawn as dark blue dots. The light blue dots 
show the corresponding radial cumulative number distribution of M dwarfs with emission. 
The solid lines with slope 3 indicate the expected radial cumulative number distribution in 
a homogeneous sample of stars. 
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Fig. 3. — The ratio of the velocity weighted number of M dwarfs with emission divided 
by the velocity weighted number of all known M dwarfs at a given age of the Galactic disk 
after its formation. We assume a present age of the Galactic disk of T = 10 Gyrs. This is a 
measure of the difference of the surface density of the disk today minus the surface density 
at time t, which illustrates the growth of the stellar disk of the Milky Way. The error bars 
are explained in the text. The solid line is a fit of relation (3) with an index of (3 = 2.2 and 
a coefficient a = 18.7. The thin dashed line illustrates the growth of the surface density at 
a constant star formation rate. 
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Fig. 4. — The correlation of the star formation r ate with th e gas s urface density. The symbols 



are data of 61 normal galaxies reproduced from iKennicutt 



(1l998l ): Sa - Sb galaxies are coded 



as open triangles, Sb - Sc as open circles, and Sc - Sd galaxies as filled circles, respectively. 
The position of the local Milky Way determined in this work is indicated by a red star. 
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Fig. 5. — The relative star formation history of the Galactic disk as determ i ned i n this work 



(circles, same as in Fig. [3l) compared with det e rmina tions by 



de la Fuente Marcos fc de la Fuente Marcos 



(rhombuses) , iRocha-Pinto et al 



3 



Cignoni et al 



(120041 ) (pentagons) 



(I2OOOI ) (downward triangles) , and 



fl2006h (squares) 



Hernandez et al. 



Vergely et al 



( 200 



(I2OO2I) (up- 



ward triangles). 



